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I. INTRODUCTION 

In the inflationary scenario graviton (metric tensor) and density (metric scalar) pertur- 
bations arise from quantum fluctuations |l[] early in the inflationary era. In many models the 
subsequent development of these separate perturbations is governed by the same parameters 
of the given model. Thus the absolute magnitude of the perturbations presently observed in 
the CMBR and through other phenomena is predicted in terms of these parameters so that, 
among other things, a comparison of the metric tensor and metric scalar contributions can 
be made. Models with power law inflation are of this type. It is the purpose of this paper 
to evaluate and compare the two contributions for such models. In doing so we shall make 
the instantaneous approximation for the transition from the inflationary era to the radiation 
era, since it holds with good accuracy over the range of metric perturbation wavelengths of 
most physical interest [0,Q ; we do not make that approximation for the transition from the 
radiation era to the matter era. 

We take power law inflation to be driven by a single canonically normalized scalar field 
with Einstein gravity. The assumption of a single scalar field is more widely encompass- 
ing than it seems; almost any current extended gravity theory, such as a higher order or 
scalar-tensor theory can be rewritten as Einstein gravity using a conformal transformation. 
Then, moreover, slow roll inflation solutions (giving the standard calculations of density 
perturbations) can be expanded about power law inflation solutions [|]. 

We use a synchronous gauge for the formalism of the computing programs, while ex- 
tracting gauge invariant quantities 0,0] for making physical comparisons and for delicate 
procedures of conveying information from one cosmic era to another. Thus we eliminate any 
confusion from unphysical synchronous gauge modes from our results. 

II. THE INFLATIONARY ERA 

We calculate with perturbations in a synchronous gauge so that the metric is 

ds^ = —a^dr'^ + a^{5ij + hij)dx^dx-' (1) 

In the inflationary era , with scalar field (p, the energy-momentum tensor is 

T,. = 0,M0,. - 9,. [^^"^0,a0,/3 + V{<P)] (2) 



A. Metric Scalar Perturbations 

In this section we shall first treat the basic formalism and then specialise to the case of 
power-law inflation. 

We adopt the notations of Grishchuk |^ so that for density perturbations with wave 
number k 

hij = h{T)QSij + hi{T)k''^Q^ij (3) 



where Q = exp(zk.x). 

If pi and pi are the perturbations to the energy and pressure densities then the Einstein 
equations yield, where n = SnG, 

a?Kpi = 3ah' + k'^h — ah'i, (4) 

a^npi = —h" — 2ah', (5) 

and also with the hypothesis of non-diagonal space-space components in the energy- 
momentum tensor, as with a scalar field or with a perfect fluid 

= h'l + 2ah'i - k^h (6) 

where a = a' /a. 

These equations being valid with an energy-momentum tensor of a scalar field or a perfect 
fluid (or a mixture) are applicable through all the model eras we shall consider from inflation 
to the present. 

We now consider the inflationary era case with the energy-momentum tensor given by 
Eq.(^. Then if the unperturbed value of is 0o and the perturbation to it is 0i the Einstein 
equations yield 

/.01 = h'/<P', (7) 

and also p 

/i" + /i[P - {a^)"/{a^)] = (8) 

f^ = ^{h' + ajh), (9) 



7 = 1- a'/al (10) 

For the particular case of power law inflation the scale factor 

a{r) oc (r, - rf (11) 

where p{< — l),rj are constants so that 7 = (1 +p)/p = constant. Our model is one of 
distinctly power-law, rather than exponential, inflation; that is we do not approach the de 
Sitter limit, p ^ —1 and 7—^0. If t is the appropriate cosmic time for the inflation era 
then a{t) oc f where q = I/7. An adequate amount of inflation requires q > 10; for larger 
q we limit ourselves to values where expressions such as Eq.(^) below are extremely good 
approximations and this is certainly true up to g = 100. 

The solution for h is given in terms of that for /i by Eq. (^ as 

a 1"^ 

h = -{Vl /i(fc,r')rfr' + Q}, (12) 



where Ci is an integration constant subsuming, when 7 is constant, the lower hmit of the 
integration, h = -Ci is a solution of /i = corresponding to an unphysical mode of the 
synchronous gauge and we drop this term since, giving zero contribution to the gauge 
invariant variables, it does not carry information through into the radiation era , as related 
in the next section. 

All the above development was non-quantum mechanical. We denote the corresponding 
quantum field theory quantities by a tilde: 

h= -v/7 / ilik,T')dT' (13) 

(Jj ^ 

— -^-==[cfc/ii(A;, r) exp(ik.x) + h.c] (14) 

zit]2\/ 2k 



{27T)2V2k 

where c^ is a quantum annihilation operator, [c^, c^/] = (5^(k — k'), and fii{y),y = |A;(rj — r)|, 
is that solution of the Bessel equation (§) (with a, 7 specified by power-law inflation as above) 
such that 

I2i{y)^e-'y,y -^00 (15) 

thus corresponding for large k to the usual mode function of quantum mechanical plane 
waves. 

A^ is a normalization factor whose determination gives the absolute magnitude of the 
observed density perturbations (in terms of the model parameters) given the assumption 
that these come from a primordial vacuum with zero quantum occupation number. The 
determination of N is as follows. 

Using the methods of ref. and Eqs.(^ ^) the gauge invariant scalar field, 0f* = 
(pi — ^(j)Qk~'^h'i is given by the gauge invariant version of Eq.(|^ as 

^^<f>T = -y7[-/i + a-fk-'{-fi' + fi{a + 7727)}]/0[, (16) 



where, for power law inflation, v^^o ~ ~Oiy/2rf. So in the limit A; ^ cxd it follows that 
-\/2^0f* — s> jl/a = Nfli/a. This yields the result 



N = y^ = VlGiiG (17) 



since then 



f ^ a W - — 3 [ckfii{k,T)exp(ik.x.) + h.c] (18) 



d^k 

{27r)IV2k 

which is the appropriate quantization for a scalar field with the normalization (§). This is 
the same normalizing factor as that found by Grishchuk using a similar argument. Also 
it is the same normalizing factor as that appropriate for the metric tensor case |0,^ given 
below. 

We can now proceed with the evaluation of the metric scalar components at the end of 
the power-law inflation. The exact expression for the solution fii is, with n = ^ — p, 
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My) = Y Y*^"^" ~ ^"^"^ exp[-i{-mi + -tt)] (19) 

As previously stated we shall assume a sudden transition at the end of inflation and the 
beginning of the radiation era; variables at that interface we denote hj t = T2, a = a2, H = 
H2 = —p{{Ti — T2)a2)~^ = (T2a2)~^. The values of k of interest are of the order of magnitude 
(the relevant meaning here is being within a few factors of ten) of 

ki = a[/ai = aiHi (20) 

where ai,Hi denote the values of the scale factor and Hubble at the time, ri, when the 
matter era begins. Thus for such values of k the values oi y = \k{ri — T2)\ at the end of the 
inflationary era, beginning of the radiation era, are of order 10~",n > 10. So then we can 
expand yU in a power series with leading terms 

..(.) = M(p)[y'-5j|^ !■ (211 

M{j,) = -2-T(l - p) exp(-j(l - p)/2)/v^. (22) 

The corresponding expansion of h from Eq.(|12|), with Cj = 0, is 



h = (J^a)-^M(p)[yP - , ^^ C^ 1 (23) 

We shall need these expressions at the interface of the inflationary and radiation eras, r = T2. 



B. Metric Tensor Perturbations 

The mode function ^/a for gravitons is given by 

/ + /i(fc2 - a" /a) = 0. (24) 

This holds for any cosmic era, whatever may be the dynamics responsible for the particular 
form of a{T). The quantum field theory expression for tensor perturbations is 0,^ 



2 ^3^ 

hij = VlQirG J2 / ,^ ,3 /— , A a\keij{k)fi{k, r) exp(zk.x) + h.c], (25) 



d^k 
t^J (27r) iv^a(r) 

where a\k is the annihilation operator for the graviton with polarization A and wave number 
k, 

/i = /ii (26) 

and the polarization tensor satisfies Y.i,j^ij{k)e^j{k) = 25\\i. 



III. INTO AND THROUGH THE RADIATION ERA 

We take the instantaneous approximation to the transition from inflation to radiation. 
This is well justified because for the values of k of interest (discussed above in [II A| ) and with 



a time of transition of order Tj — r2 = p/a2-f^2; then the parameter measuring the suddenness 
of the transition is, from Eg. (pOD , of order ki{Ti — T2) = paiHi/a2H2 ~ O2/01, and is thus 
exceedingly small. We use the method of Deruelle and Mukhanov to match the end of 
inflation to the beginning of radiation. 

Our convention is that the conformal time r is continuous from inflation through the 
matter era, and that in the pure radiation era the scale factor 

a (XT (27) 

It follows that in the radiation era, r = {aH)~^. In the inflation era the scale factor is 
proportional to (tj — r)^ (Eq.|Tl|) and in the pure matter era to (r — r^)^ where Ti,Tm 
and the constants of proportionality are determined from the continuity of a, a', with the 
convention a(present) = 1. 

A. Metric Scalar Perturbations 

In the radiation era, treated as a relativistic perfect fluid phase 

zy" + ipz/ = 0,z/ = -(/i' + «7/i), (28) 

3 a 

where 

h=- r udr, 7 = 1- a'/a'^ = 2. (29) 



a 

The mode function /i of the inflationary era is replaced by u. We write the solution of 
Eq.dD as 

u = B+ cos(fc(r - T2)/V3) - B_ sin(fc(r - T2)/V^). (30) 

The gauge invariant function $ |^ is given by 

$ = /i - {a/k^)h[, (31) 

and using the Einstein equations, Eqs.(^J^), with pi = pi/3, and Eq.([28|) 

^ = -{u'/a-u)/ae^,e = k/{aV3). (32) 

The matching conditions found in ref. 0, by deduction from the Lichnerowicz conditions 
IP, imply that $ and the expression 

r = {<^'/a + $ + e2$)/7, (33) 



should both be continuous on the surface t = T2. On the radiation side of the surface, 

T = r2+: 

<^={e2B_+B+)/{a2el), (34) 

r = -[5+(l - el) + e^B^l - \el)]/ia,el). (35) 

We now need the values of $ and F on the inflation side of the surface, r = r2_. In the 
inflationary era the Einstein equations yield 

h[ = [h" + h\a - 2a' /a - i h) + k^h]/a (36) 

Using Eq.(^) we find that to leading order at r = T2- 

$ = ^^S, r = -^—T., (37) 



S = M{p)yd^/a, (38) 

where the expansion parameter 1/2 is given by 

1/2 = kCTi - T2) = -pk/a2 = -pkT2 (39) 

noting that the continuity of a and a = a' /a is part of the matching conditions [0] . 

Solving the simultaneous equations got by equating the different expressions for the <l>, F 
pair at T2- and T2^ gives 

B+ = -62^- = 2a2S, (40) 



B- = 2p^^M{p)yr\ (41) 

and the extreme smallness of €2 means that we can put B^ = and thus neglect the cosine 
term in Eq.(|30|) for z/. Thus Eqs.(^,^) determine u and in consequence the gauge invariant 



amplitude $, Eq.(^2|), throughout the radiation era: 

$ = 5_[ecos(A;(r - T2)/V3) - sin(A;(r - T2)/V3)]/{ae^). (42) 



B. Metric Tensor Perturbations 



For gravitons the matching is simpler |]10|,3. All that is required is that hij and its first 
time derivative be continuous. 

We have, in both the inflation and radiation eras, Eq.(]25|) for h^j with /i = /ii for inflation, 
but in the radiation era given by fi" + /i(A;^ — a" /a) = fi" + k'^fi = so that 

fi = G+ cos{k{T -T2))- G. sm{k{T - ra)). (43) 

Matching (fi/a) and (/u/a)' we find, to leading order in 1/2, 

G^ = 0,G_=pMip)y^,-\ (44) 

For power law inflation p = —{1 + 6) where < 5 < 1/10. We can make a direct comparison 
of Eq.(PD with Eq. (^T]) , which would suggest that the influence of density perturbations on 
the CMB will be greater than that of gravity waves. However there are more stages to go 
through. 

IV. THE RADIATION TO MATTER TRANSITION 

We shall assume that for r > ri the universe can be characterized by matter with zero 
pressure for some value ri of the conformal time. The transition time, Ttr, from the radiation 
era, characterized by p = |p, can be as large as ~ (aiHi)^^. Thus for waves influential in 
the CMBR, section ( |1IA|) , krtr is not necessarily small and then the instantaneous transition 
approximation is not reliable. For the development of the matter components there have 
been many detailed studies [] using the Boltzmann equations for various components of the 
matter which might be present. These details are important for the smaller scale structure 
of the CMBR, but not for the larger scale structure. We wish to concentrate on this larger 
scale in our comparison of the gravity wave and density perturbation contributions, and 
thus to avoid questions of the components of matter. Consequently we shall approximate 
the transition from radiation, relativistic matter, to the non-relativistic matter domination 
era by a smooth change in an overall equation of state. The most significant parameter of 
such an approximation is the length of time that the transition takes, and while we can 
indeed make a reasonable estimate of this time, we shall also consider results as a function 
of the transition time. 

A. Metric Scalar Perturbations 

We postulate a smooth transition from the radiation era to the matter era in which we 
first parametrize the density as a function of the scale factor, a, by 

P = Qie-'\ (45) 



^See White et al. |11|, Ma and Bertschinger |jT^ and references therein. 



where s = s{r),r = \n{a/ai),ai = a{Ti), qi = constant. The energy conservation equation 
gives the pressure and thus an imphcit equation of state by 

^{pa') = -3pa^ (46) 

P/P=3(^^ + s)-l- (47) 

s = 4 corresponds to the radiation era, s = 3 to the matter dominated era; we can postulate 
an exphcit form for s as a sufficiently smooth and smoothly joining function of r between 
those two constants. The one we adopt is given in the appendix. Given such a function we 
can find a = r' = a' /a and the Hubble parameter H = r'/a as functions of a, and also, 
consequently, r = T{a) or equivalently a = a^r), when given also the values of if = Hi, a = 
ai at the beginning of the pure matter era. Thus from the above equations: 

gi = 3iHiaif/K, (48) 

a = r' = Hiai exp[(2 - s)r/2]. (49) 

Adopting the notation that the transition begins at a = ae,r = rf,,T = Xg, we specify the 
length of the transition by 

rtrans = ^{ri - Tg) = - In^ai/tte). (50) 

We now consider the density perturbations; these are given by the development of z/, defined 
in Eq.(^), or equivalently of u where 

a dh 

u = —u = a{- — l-'jh). (51) 

a dr 

In the radiation era v took a simple sinusoidal form but now evolves according to 

d^u du,^ ^ dB. .,, , ,0 
-— + — 3 + C+— + u[ikcs/ay 
dr'' dr dr 

+ (l + -)(C+l) + 2^ + 2(-n=0. (52) 

where c^ = Pq/pq,B = (2 — s)r/2,C = Scj. — -r^/c^, with po and po, the unperturbed 



du 
' dr 



pressure and density, given by Eqs.( |45|j47| ). The initial conditions for solution are that u 
are continuous with the corresponding radiation era quantities; the equation has to be solved 
by numerical methods and the result is that the physical information on the development 
of the perturbations is passed continuously from the radiation era to the pure matter era. 
From the evolution of u by Eq. (|52D , the evolution throughout the transition of the gauge 
invariant metric scalar perturbation, $, can be computed using Eqs.( 
In the matter era a oarf where 
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^ = h-^,h\ = lC,--^^C,^ (56) 



r] = T — Tm,Tm = constaut, (53) 

with the pressure being equal to zero. Consequently Eqs.(|^, ^ yield 

h = C, + C^r\ (54) 

h\ = Cieri/h + C2vlh'' + Cmk''v~V'2, (55) 

where Ci,C2,Cm are constants; the synchronous gauge mode, proportional to Cm, being 
non-physical it is just Ci, C2 that we require. The gauge invariant perturbation, $, is given 
by 

^h' = ^Ci - M 

k"^ 5 k"^!]^ 

and we note that Cm does not appear. At the beginning of the pure matter era (r = ri, r^ = 
r]i,a = ai) 

^n) = ^C^-^ia,/kYC2, (57) 

$'(ri) = ^a^{a^/kyC2 (58) 

and the continuity of $, $', calculated through the radiation and transition eras, gives Ci, C2 
in terms of (i) the inflationary power p and T2 (or equivalently 02) by Eqs. (|41| , |42|) and of (ii) 
the parameters of the radiation to matter transition. 

We can also investigate what the sudden transition approximation gives. Analogously to 



section [III A| we then have to enforce the continuity of $, F at r = ri; F is given by Eq.fpSD 
where now e = k/aVS is evaluated at r = ri, e = ei. From this we find that if we denote 
Fi_, <l>i_ to be the values at the radiation side of the interface then 

Ci = $i-(l-^e?) + Fi_ (59) 

C2 = ye?{-^<|.i_(l + e?) + Fi4 (60) 

Thus for the coefficient of the growing component we find 

Ci = {B^{2cosy — ei siny) — B^{2smy + ei cosy)} /6a (61) 

where y = kiji — r2))/-\/3 ~ krxj \pi = ei. Here B+ = is, as we have found above in 



Eq. (pOD , a result good to many powers of 10 relative to -B_. So we write 

Ci = 5_(2siny + eicosy)/6a (62) 

and in the discussion section we shall compare this with the corresponding result of ref . . 
An essential difference between this radiation era to matter era transition and the in- 
fiation to radiation era transition is that now we cannot necessarily expect, for values of k 
relevant to the CMBR fiuctuations, that the sudden approximation is nearly exact. We inves- 
tigate this quantitatively below, taking account of both the growing, Ci, and the decaying, C2, 
mode. 
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B. Metric Tensor Perturbations 

Throughout all cosmic eras in the FRW universe, gravitons of primordial origin are given 
by Eqs. (p^j25|) . Thus their development depends only on the evolution of the scale factor of 
the universe through the function a" /a. From the beginning of radiation to the present time 
the gravitational mode function /i, or equivalently the Bogoliubov coefficients p|,pr3|, evolve 



smoothly by Eq.(^J) given any twice differentiable scale parameter a{T). In the radiation 
era we have Eq.(^3|) for /i while in the matter era where a" /a = 2?7~^ we express the solution 
in terms of spherical Bessel functions: 

yi=^j—[G,J-^{z) + G2J^i{z)l (63) 

where z = krj and Gi, G^ are constants. In the radiation to matter transition // is found by 
computation using the postulated smooth aij) (see Appendix) and thus Gi,G2 are found 
in terms of the G. of Eq.(H3) through the continuity of /i, /x'. 



V. THE SACHS- WOLFE EFFECT 

Both the scalar and tensor perturbations give rise to perturbations in the wavelength of 
the photons of the CMBR. The observation of apparent temperature fluctuations, to which 
the perturbations would give rise, is well established and continuing. Scalar perturbations 
are usually considered to be dominant and we shall first deal with these. 

To be appropriate for observations nearly in the rest frame of the earth the calculation 
should be done in a comoving frame implying, 0, Tq = —a^'^h'Q^i = 0, and from Eq.(^4D 
this requires Gm = in the matter era. A synchronous coordinate system which is not 
comoving can be changed to comoving while remaining synchronous. So we thus complete 
the definition of our coordinate system, which is continuous from the beginning of radiation 
to the present time. This makes no difference to our results for Ci, G2 in terms of the -B+, -B_ 
determined at the beginning of the radiation era, as can be seen from the discussion at the 
end of |1V A|. In the matter era h and hi are now given by Eqs.(El,B^) with Gm = 0. 



We consider reception of the CMBR at the present time,?] = ?7o, and emission at 77 = 77^ 
where tje > r/i is a time within the matter era. (For the matter era we use the more 
convenient conformal time of Eq. (|53|) ) . Then in the synchronous comoving coordinate system 
the CMBR fractional temperature variation as a function of the direction of observation 
specified by the unit vector e* is |jT^ in quantum mechanical form 



AT 1 r(^E dh- ■ ■ 

T^ ' 2 Jo dr] 

I G f^E f d^k , dhik.e , .. . .n , /^,v 

= "V^io 7 7iI^^'=rf^l^^^P(^k.x) + c.c.}]rf.. (64) 

In the above equation, uj = rjo — rj, the integration is along the light path x* = e*a; and the 
quantum mechanical annihilation and creation operators, Ck and c|,, are the same as those 
in Eq.(|l^). Only hi appears in the mode function in Eq. (|64D because dh/di] = 0. 
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The angular correlation function for different directions ei and 62 is 



K={0 



6T 6T , 
^(eij — (62, 



0), 



(65) 



where |0) is the vacuum state, appropriate to the operators Ck, of the universe at the begin- 
ning of inflation. We are working in the Heisenberg picture where the cosmic state vector 
is constant and the operators vary; in Eq.(|^) the variation of the quantum operators is 
subsumed in the variation of the mode functions. 
Grishchuk [^ has shown that, with cos 6 = 61.62, 



K = K{cos6) = J2KiPi{cos6), 
1=0 



(66) 






(67) 



{( 



/(/-I) 



x^ 



-1)X '2Ji+i/2{x)+2x 2Ji+3/2(x)}. 



(68) 



We can express this with more information given as follows: We have shown that the 
parameters Ci, C2 of hi are linear functions of 5_|_, B_, say 



Ci = Ci_-B_ + Ci+-B_|_; C2 = C2--B- + C2+-B+. 

That is, for example, Ci_ is the value of Ci when calculated for B_ 
from Eqs. (^ , ^lD the value of -B_, B^ being nearly zero, so that 



1,^4 



(69) 

0. We know 



Ki = (2/ + l){^^f{^fni\p)^ r ^A-kr,pr^ 
T2 Ti p+ 1 Jo 2k 



X 



P 

'"dx[\c,.z' + C2.{^)']-Li{x) 

00 z z 



(70) 



where z = ki], x = krjQ — kt] = kuo. 

A similar procedure can be carried through for the gravitons, having the mode function 
of Eq.(|63D in the matter era with Gi = Di_G^ + D1+G+; G2 = D2-G^ + D2+G+. Knowing 
the value of G_, and that G^ is approximately zero , from Eq.(0) we find ^ 



Ki = {21 + !)/(/ + !)[/(/ + 1) - 2]{^-^n-r'M\p) r fA-knpf^ 

T2 T"! JO 2k 

I 1 2 

l-nz 1 5 



X 



rxE 

I dx\Di^Ji 



5/2{Z) - D2-J-5/2{z)]\l Y^^ 2J/+1/2(X) 



(71) 
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VI. RESULTS 

To produce numerical results we have to specify parameters giving (i) the length of the 
transition (ii) the end of the transition when the matter era begins (iii) the background 
radiation emission time for the Sachs- Wolfe calculation. We can conveniently specify using 
the scale factor, a, or the red shift, 2, where z + 1 = ao/a; our convention is that oq = 
a (present) = 1 

For (iii) it is rather natural to take the time of last scattering to be at the beginning of 
the matter era and this we shall do. We can make a remark from the literature supporting a 
not too great sensitivity to this assumption. If the emission time were taken to be within the 
radiation to matter transition then the correction from including the integrated Sachs- Wolfe 
effect is expected to be about 4% in the quadrupole moment from scalar perturbations and 



less for other moments 16 



For the begining of the matter era (zero pressure) we consider two values, mainly zi = 10^ 
but also with some illustrations from zi = lO'^'^ to assess the sensitivity of the results to the 
choice of zi. As stated above Zi is also taken to be the emmision red-shift in the Sachs- Wolfe 
calculation: 

ZE = zi (72) 

The parameter for the length of the transition, from a = Oe to ai, is given by r trans of Eq.(pUD 
so that rtrans = \ln{ai/ae) = \ln{ze/zi). We shall present some results for a range of values 
of Ttrans including Ttrans = 0, that being the sudden transition approximation which has 
often been used. We consider that a likely physical value is round about 5 e-foldings: 

f trans ^ ^-0 l'"-'j 

This can be very roughly estimated from a likely range for zeq-, which can then be used to 
make an estimate from the known [^ analytic two fluid model for the transition. 

Our results of course also depend on p, the power of inflation , Eq.(|lT]). We present 
results for a range of values of p, and firstly for the ratio of the gravitational to the density 
perturbation contributions to the multipoles. We emphasize that the ratio is calculated with 
the same physical model for both contributions and with no significant approximations. For 
the ratio no other parameters enter other than the ones just stated. But in the absolute 
magnitude of each contribution, Eqs. (|70| , [7TD , there is a common factor 



F2 = i-n-r' (74) 

containing the parameter T2, the time when the radiation era begins. The choice of T2 finally 
specifies the absolute magnitude of the multipoles. Within the model p, T2 are the important 
and significant unknown parameters; rtrans and Zi also have importance, as will be shown, 
but more is known about the possible values of these. 

As outlined in |T| our calculation of the radiation to matter transition is good for larger 
angular scale. There is an approximate relation between the correlation angle 6 and the 
order of the multipole, /, which would give the main contribution to such a correlation \W\: 

13 



We give results for multipoles with / < 8 corresponding by Eq.(^) to angular scales 9 > 8°. 
First we present the ratios, {T/S)i, of the gravitational wave to the density wave Legendre 
series coefficients, Ki, of Eqs.(^,^); these ratios are independent of T2. In Table 1 we show 
these for 4 different values of p with rtrans = 2.5 and ze = Zi = 10'^. 

From Table 2 we see,that for any particular multipole, the variation of T/S with p is 
approximately {p + l)/p. This is not quite unexpected as the ratio of gravitational to density 



wave amplitudes in the radiation era is given by Eqs. (|4l|j4^ ) as J{p + l)/12p. However Table 
2 shows that the ratio in observable multipoles, {T/S)i, is bigger than {p+l)/12p by a factor 
of order 100. 

Production of an adequate amount of inflation requires the cosmic time power q > 10. 
From Table 1 (and taking account of the proportionality to (p + l)/p = 1/q) we see that 
for q < 100 gravitational waves make a non-negligable contribution and that for q of order 
10 the contribution is similar to that of the density perturbations. Gravitational waves are 
proportionally most important in the quadrupole and octopole moments and Table 3 shows 
that these are the largest multipoles. 

Table 3 exhibits the absolute magnitude of the density component of the multipoles, from 

which the graviton induced component can be found using Table 1 . Since Ti is approximately 

known the most important adjustable parameter of the model here is T2, the conformal time 

at the beggining of the radiation era ~ end of inflation. In the limit as p — *> — 1, F2 — > (-^^). 

'''2 

We can find an order of magnitude of T2, or equivalently the more convention- free red- 
shift at the beginning of the radiation era, Z2 + 1 = ao/a2, by taking the estimate from 



observation of the quadrupole moment [jrT],|T7|. If we insert Qrms-ps ~ 18 jjK then with 



p = —1.112, zi = 10 we find, where h is the Hubble parameter 

Z2 ^ 2.2 X lO^^/i-^^ (76) 

and inserting Q/iK gives 1.3 instead of 2.2 in the above |. 

Table 4 shows the absolute values in units of 10^F2 of the scalar and tensor quadrupoles 
for a range of values of the transition length rtrans and p = —1.112, q = 10. The values for 
f^trans = wcrc obtained using the same sudden transition methods as those for infiation to 
radiation in |TTT|. For rtrans > 0.4 we use the smooth model transition of the Appendix. (For 
< rtrans < 0.4 the results are model dependent.) We note the approximate constancy of 
these results when the number of e-folds (= 2rtrans) is between 1 and 6 but that the sudden 
transition results are somewhat different, T/S then being about 3/2 larger. If we were to let 
the range of values of k/a to be even smaller than those relevant for the CMBR observations 
we would expect these numbers to approach nearer to equality. Using the comoving gauge 
Lyth [|1^] showed rather generally that, for small enough values of k/a, the development of 
the density perturbation from one well-defined cosmic era to another is independent of the 
type of transition provided it be reasonably well behaved. 



^The amplitude Q is related to the correlation multipole hy Q = Tq^/K^ [pi]]. 

14 



All the results discussed above were for values of the relevant red-shifts given by ze = 
Zx = 10^. If we increase this value to 10^'^ the scalar and tensor contributions to the 
multipole moments increase by a factor of around 4; an increase is expected as the Sachs- 
Wolfe calculation is applied to a longer photon path. The ratio (T/S)2, with p = —1.122, 
changes from 1.3 to 1.2 with similar changes for the other multipoles. 

VII. DISCUSSION 

In this work we have used the same power law inflationary model to evaluate and compare 
the density perturbation and gravitational wave contributions to the CMBR fluctuations; the 
transition from radiation domination to matter domination has been appropriately treated 
as a gradual rather than a sudden transition, but avoiding particle-spectrum dependent 
details. This latter limits the precise validity of our calculation to multipoles influencing 
correlation angles greater than about 5 degrees. For power law inflation, f^, the larger part 
of the quadrupole moment comes from gravitational waves ii q < 13.5; the corresponding 
equality value for other multipoles is ~ 10. At g = 20 gravitational waves contribute about 
1/3 to correlation Legendre coefficients of / < 8 and 10% at g ~ 100. 

We followed the paper of Grishchuk ||^ by evolving using the metric perturbations, and 
also by making computations in the synchronous gauge, but our usage of this gauge, and 
our tracking of the perturbations from one cosmic era to the next, quite differ from his. In 
particular we made no significant use of non-physical synchronous gauge modes; we used 
gauge invariant variables in critical calculations of transitions and also, where a transition 
was treated as sudden, we used the Lichnerowicz matching conditions as interpreted by 
Deruelle and Mukhanov 0. 

The message delivered by our work on the relative importance of gravitational waves 
appears to differ from that of Grishchuk 0. Since the question can be raised as to where 
this difference comes from it may be of interest to identify the detailed reasons: (i) Firstly 
in ref. a special form for a continuous transition from inflation to radiation is used. In 
principle, with an appropriate choice of form, there should be nothing wrong with such 
a treatment; there is no interface and all quantities are continuous, this being indeed the 
principle we have used in the radiation era to matter era transition. Trouble seems to arise 
in ref. when the transition is taken to the zero transition time limit. This is particularly 
evident in the flnal treatment of 7 = 1 — a'/a^ which increases from the small value (p + l)/p 
(in our notation) to 2 during the transition. The nature of the final expressions obtained 
in ref. [^ for the equivalent of 5+ and B_ and other quantities make it necessary, for a 
reasonable result, to rather arbitrarily put 7 = 2 in these expressions Q The result we find 
is, on evaluation of Grishchuk's expression for the dominating constant which we call S_, 



that it is smaller by a factor J{p+ l)/p than our expression given by Eqs.(^8|),(^0D, (^Tj). 
This makes a factor of (p + l)/p smaller for the density perturbation contribution to the final 
result, (ii) Secondly there is the question of the evaluation of the constants Ci and C2 of 



^This non-rigorous treatment of 7 has also been criticised by Deruelle and Mukhanov pi. 
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the matter era, which involves the transition from the radiation era. In expression eq.(82) of 
ref. @ for Ci, this expression has been treated as sudden (analytic expressions such as those 
of ref. 1^ cannot otherwise be obtained); there are also other approximations equivalent to 
taking only the first term in eq. (|62D . Eq.(82) [H is greater by a factor 3/2 than our first 
term in addition to having the smaller value of -B_ just related above. In our opinion this 
factor 3/2 arises from the non-trivial usage of the radiation gauge mode which is determined 
in ref. [Q by an extra continuity condition not arising from Lichnerowicz matching. Thus 
we believe there are two parts of the density perturbation calculation of ref. which are 
incorrect. 

A number of papers [|1^] by various authors found results on the relative importance of 
density and gravity wave perturbations from inflation in the CMBR fluctuations. These used 
particular properties of the various models in fairly complicated careful deductive processes. 
(For a recent account of a large range of density perturbation calculations, also referring 
to gravity wave contributions to the CMBR, and with complete references, see ref. [^.) 
By contrast we have used a direct calculation following in a unitary way both density and 
gravitational wave contributions from birth to observation, in a certain type of inflationary 
model. And our conclusions are not in great disagreement with the general concensus, for 
power law inflation, of those previous papers. 
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APPENDIX: 

We adopt the following form for the function s{r),r = ln(a/ai), of Section p.V A| where 



s(re)=4,s(ri)=3, (Al) 

and s(r) is specified between r = r^ where the radiation era ends and r = ri where the matter 
era begins. In the radiation and matter eras po ex: exp(— 4r) and exp(— 3r) respectively so 
the join with the transition era can be made arbitrarily smooth by specifying that 

^ = 0;r = re,ri, (A2) 

up to the necessary m. We take m = 3 and implement s{r) by the polynomial form 

3 1 

s(r) = So + c{r — ro)[l — x^ H — x"^ x^], 

5 7 

ro = {re + ri)/2, Sq = {s^ + Si)/2 = 7/2, 
c = 35(si - Se)/16(ri - re), 
x = (r-ro)/(ri -ro). (A3) 

The degree of smoothness specified by ^^ with m = 3, ensures the smooth joining at r = rg 
and ri of the coefficients. 
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TABLES 

TABLE L The ratio (T/S)i of tensor to scalar multipoles of the correlation function, K, for 
inflation oc f^, where t is cosmic time; q = p/{p + 1). 



p + 1 


-0.05 


-0.08 


-0.112 


-0.20 


q 


21.0 


13.5 


9.9 


6.0 


{T/Sh 


0.64 


0.99 


1.33 


2.14 


iT/S)s 


0.52 


0.81 


1.09 


1.77 


{T/S), 


0.41 


0.64 


0.88 


1.45 


(T/Sh 


0.46 


0.71 


0.97 


1.59 


(T/S), 


0.44 


0.69 


0.95 


1.58 


iT/S)7 


0.39 


0.61 


0.84 


1.41 


iT/S)s 


0.43 


0.68 


0.93 


1.57 







TABLE IL ai = 


iT/S)i^, 




p + l 


-0.05 


-0.08 


-0.112 


-0.20 


Q 


21.0 


13.5 


9.9 


6.0 


a2 


13.4 


13.4 


13.2 


12.8 


as 


10.9 


10.9 


10.8 


10.6 


04 


8.6 


8.6 


8.7 


8.7 


as 


9.7 


9.6 


9.6 


9.5 


"6 


9.2 


9.3 


9.4 


9.5 


ay 


8.2 


8.2 


8.3 


8.5 


as 


9.0 


9.2 


9.2 


9.4 





TABLE III. Multipoles Si = Ki of density perturbations in units 10^-F2 


p+l 


-0.05 


-0.08 


-0.112 


-0.20 


Q 


21.0 


13.5 


9.9 


6.0 


S2 


3.41 


2.48 


2.08 


1.82 


S3 


2.29 


1.62 


1.32 


1.07 


54 


2.01 


1.40 


1.12 


0.86 


^5 


1.43 


0.98 


0.78 


0.58 


Se 


1.24 


0.84 


0.65 


0.47 


S7 


1.20 


0.80 


0.62 


0.43 


Ss 


0.95 


0.63 


0.48 


0.33 



TABLE IV. Scalar and tensor quadrupoles (units 10^ F2), of the correlation function K, and 
their ratio as functions of the radiation to matter transition length 



ffrans 


^2 


T2 


(T/Sh 





1.72 


3.44 


2.00 


0.5 


2.08 


2.78 


1.34 


2.0 


2.12 


2.81 


1.33 


2.5 


2.08 


2.77 


1.33 


3.0 


2.09 


2.73 


1.31 
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